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Goldstone and Higgs modes have been detected in various condensed matter, cold atom and particle physics 
experiments. Here, we demonstrate that the two modes can also be observed in optical systems with only a 
few (artificial) atoms inside a cavity. We establish this connection by studying the f/(l)/^ Dicke model 
where N qubits (atoms) coupled to a single photon mode. We determine the Goldstone and Higgs modes 
inside the super-radiant phase and their corresponding spectral weights by performing both \IJ — 21 N 
expansion and exact diagonalization (ED) study at a finite N. We find nearly perfect agreements between the 
results achieved by the two approaches when Ngets down even to N — 2. The quantum finite size effects at a 
few qubits make the two modes quite robust against an effectively small counterrotating wave term. We 
present a few schemes to reduce the critical coupling strength, so the two modes can be observed in several 
current available experimental systems by just conventional optical measurements. 

It was well known that a broken global continuous symmetry in quantum phases 1,2 leads to two associated 
collective modes: the massless Goldstone mode and a massive Anderson-Higgs amplitude mode 1,3,4 (for 
topological ordered phases, see 5,6 ). The Goldstone modes have been detected in a quantum anti-ferromagnet 7 , 
a superfluid 811 and also in cold atom systems 1217 . However, the massive Higgs mode is much more difficult to 
detect in experiments. Even so, the Higgs amplitude mode was detected by Raman scattering in superconduc- 
tors 1822 and by in-elastic neutron scattering in a quantum anti-ferromagnet 7,23 near its quantum phase transition 
to a valence bond solid 24 . Unfortunately, due to the Galilean invariance, the phase mode and amplitude mode are 
conjugate variables, the conjugate pair only leads to a Goldstone mode. There is no Higgs mode inside a super- 
fiuid 8 ~ 11,19 ~ 22 . There was a report that the Higgs amplitude mode was detected in a three dimensional superfluid of 
strongly interacting bosons in an optical lattice by Bragg spectroscopy 25 . Most recently, the Higgs amplitude mode 
and its decay rate were detected in cold atoms loaded in two dimensional optical lattice near the superfluid to Mott 
transition by slightly modulating the lattice depth within a linear response regime 26 . Notably, the experiment 
detected discrete natures of the Higgs modes in a trapped system. In a relativistic quantum field theory, it is the 
well known Higgs mechanism 3,4 which generates various mass spectrum of elementary particles. Although the 
various elementary particles have been discovered with the predicted masses, the original massive Higgs particle 
stays elusive until it was finally discovered with its mass ~ 125 GeV and width ~ 6 MeV in the recent LHC 
experiments 27,28 . 

In this report, we will present the first study of the Goldstone and Higgs modes of photons inside a cavity. The 
conventional route is to look at how "more is different" emerges, namely, study how various macroscopic 
quantum phenomena emerge as the number of particles gets "more and more" 1 . Here, we will take a dual point 
of view: study how the emergent phenomena evolve as the number of particles becomes "less and less". This dual 
approach becomes especially important in view of recent experiments of cold atoms inside an optical cavity 29 32 or 
superconducting qubits 33,34 or quantum dots 35 37 inside a microcavity involving only finite to even small number 
of particles (Fig. 1(a)). We demonstrate this dual approach by studying the 17(1) Dicke (Tavis-Cummings) model 
Eq. ( 1 ) where N cold atoms or qubits coupled to a single photon mode inside a cavity (Fig. 1 (a) ) . It was known that 
in the thermodynamic limit 38 " 44 , when the atom-photon coupling g is sufficiently large, the system undergoes a 
quantum phase transition from a normal phase to a emergent superradiant phase which breaks the global 1/(1) 
symmetry. We perform both l/J = 2/N expansion and exact diagonalization (ED) study on how the Goldstone 
mode and Higgs amplitude mode inside the superradiant phase evolves as the N decreases to a few. We find that 
even for a few number of N, the system's energy levels in the super-radiant phase display a Landau-level like 
structure with the inter-Landau energy scale setting by the Higgs energy E H and the intra-Landau energy scale 
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(a) (b) 

Figure 1 (a) AT (artificial) atoms are placed on anti-nodes of a cavity, u is the repulsive qubit-qubit interaction which can be tuned to reduce the 
critical coupling g c (b) The analytical Mandel factor Q M (red) against the ED result (blue) at N = 3. It is a number squeezed state inside the superradiant 
phase. 



setting by the Goldstone energy E G . In both the photon and photon 
number correlation functions, we evaluate the low frequency 
Goldstone mode E G , the high frequency Higgs mode E H and their 
corresponding spectral weights C G and Cu- The Higgs mode is a 
sharp mode protected by the (7(1) symmetry at any finite N. We find 
nearly perfect agreements between the results achieved from the II] 
calculations with those from the ED studies in all these physical 
quantities even when N gets down even to N = 2. We also study 
the effects of the counter rotating wave (CRW) term by the II J 
expansion and find that the finite size effects for a few qubits N ~ 
2-5 dominate those of the CRW terms if g'/g < 1/3. We discuss 
several schemes to reduce the critical coupling, so the two modes 
can be observed in several experimental systems by conventional 
optical detection methods such as the florescence spectrum measure- 
ment 45 on Eq. 5 and the HanburyBrown-Twiss (HBT) type of mea- 
surement 46 on Eq. 6 respectively. 

Results 

Reducing the U(l)/Z 2 to the / - U(l)/Z 2 Dicke model. In the t/(l)/ 
Z 2 Dicke model 47 , a single mode of photons couple to N two level 
atoms with same coupling constants g and g' . The two level atoms 
can be expressed in terms of 3 Pauli matrices a a ,a = 1, 2, 3. The 1/(1)/ 
Z 2 Dicke model can be written as: 

N N 

h u(1)/Zi =^ a + a -^zZ aZ ' + -kzZ ( fltff r + hx ) 

(1) 

+^X>V +*•*). 

where the co a , co^, are the cavity photon frequency and the energy 
difference of the two atomic levels respectively, the g = \/Ng is the 
collective photon-atom coupling (g is the individual photon-atom 
coupling). The g' = VNg' is the counter-rotating wave term. It was 
demonstrated in 48,49 that in the thermal or cold atom experiments 29,30 , 
the strengths of g and g' can be tuned separately by using circularly 
polarized pump beams in a ring cavity. In the qubit 33-34 or quantum 
dot 35 37 experiments, the CRW terms and RW terms have the same 
strength at the bare level, however, the CRW term is usually much 



smaller than the RW term at the effective level as is the case in the 
experiment 33 . This is because the former violates the energy 
conservation, while the latter respects the energy conservation, 
while the latter respects the energy conservation. However, when 
the coupling strength gets close to the the transition frequency, the 
CRW term becomes comparable to the RW term as is the case in the 
experiment in 34 . In any case, the Hamiltonian Eq. (1) with 
independent g and g' is the most general Hamiltonian describing 
various experimental systems in various coupling regimes under 
the two atomic levels and a single photon mode approximation. 

One can introduce the total "spin" of the N two level atoms 
/ z = ^ o], /+ = ^ ff+, 7" = ^2 a~, When all the N atoms 
are in the ground state, then / = N/2, J z = —N/2, because the total 
spin J 2 =Jx+Iy+ ll is a conserved quantity, by confining the Hilbert 
space only to / = N/2, then one reduces the Hilbert space from 2 N to 
2/ + 1 = N + 1. One can call the resulting model as the / — U(l)/Z 2 
Dicke model. 

One main advantage of this reduction is that one can study the / — 
U(l)/Z 2 model by using Holstein-Primakoff (HP) representation of 

the angular momentum operator J z = Vh -/,/+ = $ \J 2] — b^b, 

/_ = \J 2/ — b^bb, therefore treat photon and atom on the same foot- 
ings. This advantage will enable us to bring out many new and 
important results hard to retrieve from the 1/N expansion in 43 . 
Very fortunately, this reduction will not change the most important 
physics of the original U(l)/Z 2 model Eq. (1). As argued in the 
Methods section, except the U(l)/Z 2 Dicke model contains some 
additional energy levels, both models share the same other physical 
quantities to be studied in this article. 

lfg' = 0, the Hamiltonian Eq. (1) has the 1/(1) symmetry a — * ae' B , 
a~ —> o~e'". The CRWg' term breaks the 1/(1) to the Z 2 symmetry a 
— > —a, <r~ — » — If g' = g, it become the Z 2 Dicke model studied 
in 44 . In this article, we focus on the 1/(1) Dicke model, but will also 
consider the effects of the small counter- rotating wave term g' < g in 
the experimental detection section and the Methods section. The g' 
= g and the g' ~ g cases will be studied in 44 . The 1/(1) Dicke model 
was solved in the thermodynamic limit N = 00 by various meth- 
ods 38 ~ 44 . In the normal phase g <g c = \/03~0Jb, (a) = 0, the 1/(1) 
symmetry is respected. In the super-radiant phase g > g c , (a) 0, 
the 1/(1) symmetry is spontaneously broken. 
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Goldstone and Higgs modes in the super-radiant phase by 1/J 
expansion. In the super-radiant phase g > g c and also not too 
close to the quantum critical point (QCP) (if too close, then a'a^.j, 
b^b^ij, a direct l/j expansion is needed and will be performed 
elsewhere), it is convenient to write both the photon and atom in 

the polar coordinates a = \] ). 2 + 5p a e'°" , b = \J ll + dp b e' Bi where 

X a ~X b ~j. When performing the controlled 1// expansion, we keep 
the terms to the order of ~j, ~ 1 and ~ l/j, but ignore orders of l/j 2 . 
We first minimize the ground state energy at the order j, we found the 



saddle point values of X a and A b : X, 



COa V 2 

It holds only in the 



Xh = vW—i") where ji = g 2 c / ' g 2 = a} a u) b j 
superradiant phase g > g c . 

Observe that (1) in the superradiant phase g > g c , I 2 ~ A 2 , ~j, (2) it 
is convenient to get to the ± modes: 8 ± = (6 a ± 6 b )l2, Sp ± = <5p a ± 
dp b , k\ = a 2 + a\. (3) paying a special attention to the crucial Berry 
phase term in the 0+ sector, (4) after shifting 0- —* 0- + n/2, then 
one can get the effective action up to the order of l/j: 

£cj(i) [Sp+,0+] =i{k\ +Sp + )8 z 0 + +i(A 2 _ + Sp_)8 T 0_ 
+ ®(5p + ) 2 +D-(Sp_ +y<5p + ) 2 +4co fl / 2 sin 2 0_, 

where the first line are the crucial Berry term in the 0+ and 6- 

2(0 a g 2 

respectively, D— ^ 2 — is the phase diffusion constant, D_ = 

H CO 2 

E 2 H /\6A 2 a w a with E 2 H = (w a + (o b ) 2 +4g 2 ll/N. The y= i f(l- 

H 

s 4 

— is the coupling between the + and — sector. Under the (7(1) 
w 4 / r 6 

transformation 0 a/b — > 0 a/b + %, 0 + — » 0+ + %, 6-—>0-, so the 6 _ is 
neutral under the (7(1) transformation. There is a mass term for 0_, 
but no mass term for 0 + . The conjugate pair (0 + , &p+) leads to the 
Goldstone mode E G as shown in Eq. (5). While the conjugate pair 
(9-, Sp-) leads to the Higgs mode E H as shown in Eq. (6) (See also 
the Methods section). 

Defining the Berry phase in the + sector as a 2 + — P + a where 
P— 1,2, ■ ■ ■ is the closest integer to the k 2 + , so —1/2 < a. < 1/2. In 
fact, P = a 1 a + Vb is just the conserved total excitations number. 
Redefine 5p + =N — P, then one can write the corresponding 
Hamiltonian of Eq. (2) as: 

H, 



2 +D_[Sp_+y5p + ] 2 +4wJ 2 sin 2 6_. 



(3) 



Because the 0_ is very massive, after pinning 0_ around 6- ~ 0, one 
()_~0 2 _, so the total wavefunction is 



can approximate sin 



1 



i\(m + l)B + +y(m + l)e_] 



2n 



$i(9 J) where the / = 0,1, 



■ ■ ■ are the Landau level indices, the m = — P, — P+l, ••• are the 
magnetic indices at a given sector P, 0 < 0 + < 2n, — 00 < 0- < 00 and 
the \jji(0-) is just the Z-th the wavefunction of a harmonic oscillator. 
The corresponding eigen-energy is 

E 0 {l,m) = {l+l/2)hE H +®{m + l-a) 2 . (4) 

The ground state energy is at I = 0, m = 0. One can see that the 
energy spectrum Eq. (4) has a Landau-level structure: the Landau 
level energy scale is given by the Higgs energy E H ~ 1, the intra- 
Landau level is set up by the Goldstone energy scale E G ~ l/j. In the 
large j limit, there is a wide separation of the two energy scales 
Eh ~ 12>£g ~ 1 /j- When the excitation number P reaches the order 
of N, then the intra-Landau levels with |m| £ P will start to overlap 
with the inter-Landau levels. These analytical results explain pre- 



photon operator: a — 



+ 



At a finite N, 



cisely the ED energy level structures shown in Fig. 2 for the resonant 
case oj a = co b . 

Away from the QCP, one can write down the II] expansion of the 

SPa (Sp a , 
2l a BXl 

due to the restoration of the !7( 1 ) symmetry by the phase diffusion in 
the 0+ sector, any (7(1) non-invariant correlation functions vanish 
(a) = 0, (a(T)fl(O)) = 0 So we need only focus on the (7(1) invariant 
correlation functions. By using both canonical quantization and path 
integral approaches, we find the single photon correlation function: 

(Ta(T)« t (0)) = C G e-^ + + 0(l/j), 

C G = A 2 a -C 0 + (l-ya/2), E g = d[ 1 



C 0 = 



4E h 



-OJ b 



where E G = D 



E 0 — E H +E G , 



is the Goldstone mode with the correspond- 



ing spectral weight C G , while E a = E H + E G is the optical mode with 
the corresponding spectral weight C a . All these quantities can be 
directly measured by the florescence spectrum measurement 45 . 

The E G , C G and E m C a are compared with the ED results in Fig. 3 
and Fig. 4 respectively. One can see that except at the first few PCN 
steps, the ED in E 0 match the analytical relation E a = E H + E G in Eq. 
(5) well. The discrepancy at the first few steps is not surprising, as said 
previously, if too close to the QCP, a direct l/j expansion is needed 
and will be performed elsewhere. However, the agreement between 
the analytical and ED results in C D holds in all couplings even near the 
QCP. 

One can also compute the photon number correlation function: 

<T« a (tK(0)) - (n a ) 2 = (S Pa (z)d Pa (0)) = °^e- E "\ (6) 




Figure 2 | The ED results (See the Methods section) of the energy levels E 
measured by subtracting the ground-state energy versus glg c at resonance 
o a = (Ob with N = 5 atoms. Different colors of the energy curves 
correspond to several smallest numbers of total excitations number P = 
a f a + tfb. The dashed vertical lines correspond to the critical values of g 
where the number of total excitations P in the ground state increases by 
one. 
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(red line) are contrasted with the ED result Eg =£? 



Figure 3 | (a) The analytical Goldstone mode at y. = —1/2, Eg(j.= — 1/2) =D(g) - 

(blue lines) at JV= 5, 3, 2, 1 respectively. It is remarkable that the analytical result can even map out broad peaks at small Pin the ED results very precisely, 
(b) The analytical spectral weight (red) of the Goldstone mode Co against the ED result (blue) at JV = 3. 



where (n a ) =h a . The Higgs energy E H and the corresponding spec- 



tral weight Ch = 



E H 



are compared with the ED results in Fig. 5. 



Note that the sharpness of the Higgs mode is protected by the con- 
servation of Sp + in Eq. (3). Both Ch and E H can be directly measured 
by the HanburyBrown-Twiss (HBT) type of measurement on two 
photon correlation functions 46 . 



From the Eq. (6), one can see that \(5p a ) 



so one can 



find the Mandel Q factor: Q A 



which was compared with 



the ED result in the Fig. 1(b). For a> a = cOj,, one can see — 1 < Q M < 
— 1/2. So it is always in a number squeezed state. As 00 limit, Q M 
— » — 1 , so it approaches a photon Fock state. It is known that number 



squeezed states could be very important in quantum information 
processing and also in high-resolution and high sensitivity measure- 
ments. Very similarly, one can evaluate the atom correlation 
functions. 

Effects of the CRW term and experimental detections of the 
Goldstone and Higgs modes. The effects of the CRW terms on 
system's energy Eq. (4), photon correlation function Eq. (5) and 
the number correlation function Eq. (6) are examined in the 
Methods section. Their effects were found to be much smaller than 
those of the finite size for a few qubits N ~ 2 — 5iig'/g< 1/3. Recent 
experiments 2930 reached the Z 2 super-radiant regime 44 with the help 
of a transverse pumping. In this transverse pumping scheme, the 
CRW terms in Eq. (1) are as important as the RW ones g' = g, so 
only the Z 2 super-radiant phase can be realized. However, it was 
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(a) 

Figure 5 | (a) The analytical Higgs energy E H (red) against the ED result Eh 
for the Higgs mode against the ED result (blue) at N = 3. 

demonstrated in 48,49 that the strengths of g' and g can be tuned 
independently by using circularly polarized pump beams in a ring 
cavity. So we expect that g'/g < 1/3 can be achieved in this transverse 
pumping scheme, then the system can be tuned to the 1/(1) 
superradiant regime. It is also promising to reach the Z 2 superra- 
diant regime "simultaneously" (namely without any transverse 
pumping) with artificial atoms such as superconducting qubits 
inside micro-wave circuit cavity 33,34 and quantum dots inside a 
semi-conductor nano-cavity engraved in a photonic crystal in 
Fig. 1(a) 35 37 . Indeed, very recently, by enhancing the inductive 
coupling of a flux qubit to a transmission line resonator, a remark- 
able ultra-strong coupling with individual g~ 0.1 2iO a was realized in 
a circuit QED system 33,34 . In this simultaneous scheme, due to the 
violation of the energy conservation, the CRW term is usually much 
smaller than the R W one g' < g, but gets stronger as the coupling gets 
stronger. The effects of the CRW terms on system's energy Eq. (4), 
photon correlation function Eq. (5) and the number correlation 
function Eq. (6) are examined in the Methods section. Their effects 
were found to be much smaller than those of the finite size for a few 
qubits N ~ 2 — 5 if g'/g < 1/3. In real experiments of 
superconducting qubits or quantum dots inside a cavity in 
Fig. 1(a), there are always the potential scattering term kj^alj 
between the cavity photons and the qubits and the qubit-qubit 
interaction term u] 2 z lj. The critical coupling g c is shifted to: 

gc +g' c = \/{(a a -A z ){w h -2u), (7) 

which indicates that the two repulsive interaction terms decrease the 
critical g c well below the bare critical frequency yfw a Wj,. The qubit- 
qubit interactions can be tuned inductively or capacitively. This fact 
could be used to put the system into the regime where the CRW term 
satisfies g'/g < 1/3, so the 17(1) super-radiant phase can be realized in 
the possible future experiments using both atoms inside a optical 
cavity or qubits inside a microwave circuit QED in Fig. 1(a). 

There are also other promising experimental systems to realize the 
1/(1) super-radiant phase. Most recently, the giant dipole moments 
of intersubband transitions in quantum wells have pushed the system 
into the ultrastrong light-matter coupling regime in semiconductor 
heterostructures 50,51 . Very recent experiments 52 achieved very strong 
coupling between an ensemble of s = 1/2 spins and photons in 
electronic spin ensembles coupled to superconducting cavities. The 




(b) 

Ef — Eq (blue) atiV = 3. (b) The analytical spectral spectral weight C w (red) 



strong coupling regimes are also realized in ion Coulomb crystals in 
an optical cavity 53 . 

Discussion 

Quantum mechanics describes the motion of a single or a few part- 
icles 5458 . Condensed matter physics studies various emergent 
quantum phenomena of macroscopic number of interacting part- 
icles. Ultracold atom systems and optical cavity systems can provide 
unprecedented experimental systems to study quantum phenomena 
ranging from a few particles to a million number of interacting 
particles. Due to the tremendous tunability of all the parameters in 
these systems, they can be tuned to scale up from the isolated 
quantum mechanics systems to macroscopic condensed matter sys- 
tems. In this article, we show that the many body theory developed to 
study the emergent phenomena of condensed matter systems can 
also be a very powerful tool to study the physical phenomena from 
millions of particles down even to a few particles. Especially, we study 
how the emergent Goldstone and Higgs modes evolve as the number 
of particles gets less and less, even down only a few particles in 
quantum optical systems. The discrete natures of both modes shown 
in Fig. 3 and Fig. 5 at a finite N are due to the Berry phase effects — 1/2 
< a < 1/2. Both modes are well defined sharp quasi-particle excita- 
tions with no damping. Especially, the sharpness of the Higgs mode is 
the 17(1) symmetry protected at any finite N. We also found that the 
finite size effects at a few qubits making the 17( 1 ) super-radiant phase 
quite robust against the counter rotating wave (CRW) term. Both 
modes can be detected even with a few qubits or ions inside a QED 
microwave cavity by conventional optical measurements 45,46,59 . 
Considering it is difficult to scale up these systems to large number 
of qubits at the present technologies, this feature becomes experi- 
mentally appealing. Our theoretical works should provide a solid 
foundation for various ongoing and upcoming systems with a small 
number of particles to observe the novel phenomena due to strong 
light-matter interactions explored in this report. 

Methods 

Exact diagonalization (ED) study. For simplicity, in the following, we limit our ED 
study only to the resonant case oj a = <D b . We assume P < N. The P > N case can be 
similarly addressed by changing P + 1 to JV + 1. The ground state in the given P 
Hilbert space is: 
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|P,G) = |P,/ = 0> = ^Af f -°|N/2,5-N/2> A |P-5) f , 



(8) 



where the coefficients A^ l=0 can be determined by the ED. From Eq. (8), one can 
evaluate the Mandel Q factor Q M = — 1 + {{3n p ) 2 )l{n p ) which was compared with the 
analytical result in Fig. 1(b). 

The Z-th eigen-state in the P + 1 sector with the eigen-energy 
Ef +1 ,Z = 0,l, ---,P + lis: 



>+l,l) = ^A* +1 > l \N/2,s-N/2) A \P + l-s) F , 



(9) 



where the coefficients A^ +1,/ can be determined by the ED. 

In the Lehmann representation, we can evaluate the photon-photon correlation 
function: 



(Tfl(T)fl" r (0) ) = e~ | {P+ HW T 

l-o 

{P+l,l\a f \P,G} = A? +'X'VP+ 1 -s. 



(10) 



where Eq — Eq + 1 — Eq is the Goldstone mode with the corresponding spectral weight 
C G = \(P+ 1,1= 0\a f \P, I = 0)| 2 , while E 0 = £[ +1 -Eq is the optical mode with the 
corresponding spectral weight C 0 = \{P + 1, 1 = l^l-P, / — 0)| 2 and so on. In fact, 
there are P + 2 lines, we just focus on the two lowest energy excitations Z — 0, 1. 
Very similarly, one can evaluate the photon number correlation function: 



(Tn a (T)n a (0)} - (n a ) 2 = £ e~( E > - E ">\{P,l\n a \P,G)\\ 
i=i 

{Pj\n a \P } G) = -Y,Af' l As'\ 1>1, 



(11) 



where (n a )= f =0 1 Af '° | 2 (,P — 5) =2f ( and the Higgs mode E H = E{ - E P G with the 
spectral weight Ch — \{P,l= l\n a \P, G}\ 2 . Very similarly, one can evaluate the atom 
correlation functions. 

Relations between / — 1/(1 )IZ 2 Dicke model and the U(l)/Z 2 model. The energy 
levels in the lowest Landau level (LLL) shown in Fig. 2 are identical in (7(1) and / — 
17(1) models. This is because the ground state must be a totally symmetric state. In 
fact, every ground state in a given P — cfia + b^b sector must be a totally symmetric 
state. It is the crossings of all these ground states at different P sectors which lead to all 
the energy levels in the LLL shown in the Fig. 2. This explains why the diffusion 
constant D achieved by 1 // expansion in this article is identical to that achieved by the 
l/N expansion in Ref. 43. Because both photon and total spin operators are also totally 
symmetric in the atom operators, then all the energy levels coupled to the ground state 
by the photon and total spin operators are also totally symmetric, so this also explains 
why we achieved the same single photon or atom correlation functions in the reduced 
Hilbert space in the / — (7(1) Dicke model by 1// expansion as those in the whole 
Hilbert space by the l/N expansion. However, compared to the reduced Hilbert space 
in the / — U(l) Dicke model, there are many extra energy levels in the whole Hilbert 
space in the (7(1) Dicke model, but they are not coupled to the ground state by the 
single photon or atom operators. Similar arguments apply to the more general (7(1)/ 
Z 2 model with the CRW term in Eq. (1) and the / - U{l)/Z 2 model. 

Comparisons with the Higgs mode and pseudo- Goldstone mode in one gap and 
two gaps superconductors. It is constructive to compare the Goldstone and Higgs 
mode of the atom-photo system studied in this report with those in (charge neutral) 
superconductors (so one can ignore the Anderson-Higgs mechanism for the sake of 
explaining physical concepts). In a one gap superconductor, as explicitly 
demonstrated in the last reference in Ref. 19-22, when integrating out the fermions, 
the amplitude and phase of the paring order parameter v[/ — Ae m emerges as two 
independent degree of freedoms, instead of being conjugate to each other. Its phase 
fluctuation in 6 leads to the Goldstone mode, while its amplitude fluctuation in A 
leads to the Higgs mode. 

Now we consider the collective modes in a two gap superconductor such as M g B 2 
which has a a band and a n band. Therefore it has two order parameters ¥„■ — A a ^ 8a 
and ^-A/". There are also fermionic degree of freedoms: 0 electrons and n 
electrons. If ignoring the interband scattering V atm the Hamiltonian has two inde- 
pendent (7(1) symmetries: U(l) a X [7(1)^, the systems is just two copies of single band 
superconductor. So there are two independent Goldstone modes. 

9 a> 6 n and also two independent Higgs modes A ff , A^ for the two bands respectively. 
Now when considering the interband scattering term V^ n , the symmetry of the 
Hamiltonian reduces from U(l) a X C7(l)„to [U(l) a X (7(l)J Z) wheretheDmeansthe 
simultaneous rotation of the two order parameter phases. Then the two Goldstone 
modes couple to each other and split into one gapless Goldstone mode 9+ — 9 a + 9 n 
plus a gapped pseudo-Goldstone mode 9- — Q a — 0 n . The pseudo -Goldstone mode 
6- is just the relative phase mode between the two order parameters whose gap is 



proportional to the strength of the interband scattering V a _ n . The two Higgs modes 
A a , A n will also couple to each other and split into two new Higgs modes. In all, the 
two gaps superconductor has one gapless Goldstone mode and 3 gapped modes: one 
pseudo-Goldstone mode and two Higgs modes. 

A pseudo-Goldstone mode is always associated with an explicit symmetry breaking 
of a Hamiltonian, its gap is proportional to the strength of the explicit symmetry 
breaking. 

In contrast, a Higgs mode is the magnitude fluctuations of an order parameter. It is 
always associated with a spontaneous symmetry breaking in a ground state. The final 
physical meaning of a relative phase mode depends on the physical degree of free- 
doms of a system and its original relation to the order parameters of the system. As 
shown below Eq. (2) in the main text, the conjugate pair (t)p_, 0_) fluctuation leads 
directly to the photon amplitude fluctuation mode, namely, the Higgs mode in Eq. (6). 

To some extent, the photon-atom system studied here is similar to one gap 
superconductor discussed in 19-22 with the photon corresponding to the pairing order 
parameter, while the atoms corresponding to the fermions. When integrating out the 
atomic degree freedoms, the amplitude and phase of the photon order parameter 
emerges as two independent degree of freedoms, instead of being conjugate to each 
other. Its phase fluctuation leads to the Goldstone mode, while its amplitude fluc- 
tuation leads to the Higgs mode. This fact was demonstrated by the 1/JVexpansion in 43 
and also by Eq. (5) and Eq. (6) of this report by 1/7 expansion. As shown in the 
Methods section, a small counter-rotating wave g' term in Eq. (1) break sthe (7(1) 
symmetry to a Z 2 symmetry, then the Goldstone mode at JV — 00 will become a 
pseudo-Goldstone mode whose gap is proportional to the strength of the counter- 
rotating wave term. 

The effects of the counter-rotating wave term at N = °° and at a finite N. Now we 

consider the effects of the counter-rotating wave (CRW) terms in Eq. (1). Following 

g+g' [j 

the same procedures in the main text, we find that X a — \ - (1 — W 2 ), 

oj a V 2 

Xh — yj{} — f-t) where fi — oj a cOb/ (g + g' ) 2 , so the QCP is shifted to g+g! =g c = 
^co a cOfj. The Hamiltonian to the order of l/j is: 



Hut!)/* = -{&P+ -ot) 2 +D_ [3p_ +ydp + ] 2 



+4co a ^^-rsin 2 
g+g 



(12) 



where D— — is th e phase diffusion constant, D_ — £#/l6^0) with 



4 = K + »fc) 2 + 4(iT+s') 2 ^/W. They = j| (l- 



(S+f) 



is the coupling 



between the + and — sector. 
Eq. (12) can be rewritten as 



+2w a Xl^-- (1-OW20+), 
g + g' 



(13) 



where H U(l ) takes the same form as Eq. (3) with the parameters corrected by g' . The 
last CRW term breaks the (7(1) symmetry to Z 2 symmetry 9 a/ b — » 9 a /b + n,9 + ^> 9+ 
+ k, 9— — > 9-, so the 9- is neutral under the Z 2 transformation. In the ther- 
modynamic limit N — oo, it leads to a small mass term for 0+, so the Goldstone mode 
at N = 00 becomes a pseudo-Goldstone mode with a small gap 

A P Q — — — - — \(g +g') 4 —g^ \ ■ Obviously, this gap vanishes at the QCPg + g' — g c . 

Eh g+g' L 1 
In the following, we discuss its effects at a finite N. 

If we ignore the CRW term, all the results achieved in the main text on the systems's 
energies Eq. (4), the photon correlation function Eq. (4) and the photon number 
correlation function Eq. (6) remain intact after making the corresponding changes in 
the parameters. Then for small g'/g, at a finite N, we can can treat the CRW term by 
the perturbation theory. Here, we only list the main results. Obviously, the high 
energy Higgs mode is in-sensitive to this CRW term, so we only need to focus on its 
effect on the low energy Goldstone mode. Then the sole dimensionless small para- 
2 g' I 

meter is 3 — 2w a X a ^ f I D. (1) For the Berry phase a 0, non- degenerate per- 
turbation leads to the correction to the system's eigen-energy Eq. (4) at the second 
order — 3 2 . Note that although at a — — 1/2, the energy is doubly degenerate with 
(3p _,_ = m, 3p + — — m — 1), but m and — m — 1 carry opposite parities, so they will 
not be mixed by the CRW term. So the nondegenerate perturbation theory is valid. 
For the Berry phase a — 0, because the two degenerate states (m, — m), m > 0 carry the 
same parity, one need to use the degenerate perturbation theory to treat their splitting. 
The pair (m, — m) will split only at the m — the order degenerate perturbation, so the 
splitting AE — 3'". (2) The normal photon correlation function Eq. (5) receives a 
correction ~- 3 2 in both energy and spectral weight. Most importantly, there appears 
also an anomalous photon correlation function (Ta(z)a(O)) ~ 3 2 . So the detection of 
a small anomalous photon correlation function by phase sensitive homodyne 
experiments 45,46 ' 59 , could be used to determine the strength of the CRW term. 

One can see that the corrections to all the physical quantities are at the second order 
— 3 2 or higher. From the N = 2 qubits in the Fig. 3(a), one can see that D — coJA, 
2X 2 ~ 1 near the QCP, then when g'/g < 1/3, the corrections due to the CRW term is 
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suppressed compared to the finite size effects. Physically, at N — t», any CRW term 
will transform the gapless Goldstone mode into a pseudo-Goldstone mode whose gap 
is proportional to the strength of the CRW term. In contrast, at a finite N, the 
quantum finite size effects already opened a gap to the Goldstone mode which is of the 
phase diffusion constant D 1/N. This gap make the Goldstone in a finite system N — 
2-5 quite robust against the CRW term ifg'/g < 1/3. 
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